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ABSTRACT 

Using the idea of regularisation of singularities due to the variability of the fundamental con¬ 
stants in cosmology we study the cyclic universe models. We find two models of oscillating 
and non-singular mass density and pressure (”non-singular” bounce) regularised by varying 
gravitational constant G despite the scale factor evolution is oscillating and having sharp turn¬ 
ing points (’’singular” bounce). Both violating (big-bang) and non-violating (phantom) null 
energy condition models appear. Then, we extend this idea onto the multiverse containing 
cyclic individual universes with either growing or decreasing entropy though leaving the net 
entropy constant. In order to get an insight into the key idea, we consider the doubleverse 
with the same geometrical evolution of the two “parallel” universes with their physical evo¬ 
lution (physical coupling constants c(f) and G(f)) being different. An interesting point is that 
there is a possibility to exchange the universes at the point of maximum expansion - the fact 
which was already noticed in quantum cosmology. Similar scenario is also possible within 
the framework of Brans-Dicke theory where varying G(f) is replaced by the dynamical Brans- 
Dicke field 0(f) though these theories are slightly different. 
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1 INTRODUCTION 

One of the most fundamental problems of cosmology is the ques¬ 
tion about initial conditions which in the standard approach involve 
the discussion of singularities. The singularities appear as a result 
of our ignorance about extreme states of the universe, where one 
faces a blow-up of common physical quantities leading to the ter¬ 
mination of the evolution of the universe (at least from our point of 
view). In the standard case of ACDM cosmology a big-bang initial 
singularity appears and it marks our lack of knowledge about what 
was the state of the universe before. In general relativity, a singu¬ 
larity is defined as a region for which one cannot extend a (null 
or timelike) geodesic. This is known as geodesic incompletness 
(Hawking & Ellis 1973). However, this definition is very general 
and some more subtleties on the nature of singularities have been 
investigated leading to some milder (weak singular) states of the 
universe than the standard big-bang (strong singular) case (Tipler 
1977; Krolak 1986; Fernandez-Jambrina & Lazkoz 2004). 

It emerges that there is the whole variety of non-standard (or 
exotic) singularities (Nojiri et al. 2005; Djbrowski 2014). Among 
them a big-rip (strong) singularity associated with the phantom 
dark energy (Caldwell 2002; Djbrowski et al. 2003; Caldwell et al. 
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2003), further classified as type I, a sudden future singularity (SFS 
or type II) (Barrow et al. 1986; Barrow 2004; Freese et al. 2008; 
Shtanov & Sahni 2002), as well as numerous other types such as: 
generalized sudden future singularities (GSFS), finite scale factor 
singularities (FSF or type III), big-separation singularities (BS or 
type IV), w-singularities (type V) (D^browski & Denkiewicz 2009; 
Fernandez-Jambrina 2010), directional singularities (Fernandez- 
Jambrina 2007), little-rip singularities (Frampton et al. 2011), and 
pseudo-rip singularities (Frampton et al. 2012). All these singular¬ 
ities except big-rip and little-rip are weak. 

Singularities seem to be an obstacle to construct cyclic models 
of the universe - something which was a strong ancient Greek cul¬ 
ture believe (Phoenix Universe). In relativistic cosmology this idea 
was first developed by Tolman (1934) and further studied by nu¬ 
merous authors (Markov 1984; Linde 1984; Barrow & Djbrowski 
1995; Sahni et al. 2015). However, the crucial issue was that de¬ 
spite one was able to formally construct cyclic models, one could 
not avoid the problem of initial (strong) singularity. The finite num¬ 
ber of cycles just shifts the singularity problem into the past initial 
cycle. Besides, the physical mechanism of a transition from a cycle 
to a cycle remained unknown because of singularities joining the 
cycles. 

Recently, this problem has been approached (Barrow et al. 
2004; Alexander et al. 2016) by the application of a possibility 
that fundamental constants in the universe become dynamical fields 
(Uzan 2011). They influence the singularities in a way that these 
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become milder (weaker) or just disappear. So far mainly the theo¬ 
ries with dynamical gravitational constant G (Brans & Dicke 1961), 
the fine structure constant a (Barrow & Magueijo 1998; Barrow 
2010), the electric charge e (Bekenstein 1982), and the varying 
speed of light c (Barrow 1999; Albrecht & Magueijo 1999; Barrow 
& Magueijo 1999; Gopakumar & Vijayagovindan 2001) have been 
investigated. It has been found that varying a and varying c theories 
allow the solution of the standard cosmological problems such as 
the horizon problem, the flatness problem, and the A-problem. 

In our previous papers (Djbrowski & Marosek 2013; 
D^browski et al. 2014) we have applied varying constants theories 
to regularise (remove or weaken their nature) various types singu¬ 
larities in cosmology. This idea coincides with an idea of regular- 
isation by the higher-order curvature terms (Barrow et al. 2008; 
Houndjo 2010) which are of quantum (or stringy) origin and obvi¬ 
ously they may be equivalent to the effect of the varying constants. 
In this paper we will further regularise some other cosmological 
singularities and study the cyclic models of the universe which are 
due to the variability of the fundamental constants. We will also ex¬ 
tend this idea to the multiverse containing individual universes with 
regularised singularities which have both growing and decreasing 
entropy though leaving the total entropy constant. 

The paper is organised as follows. In Section 2 we present 
the field equations for varying c, G, and a universes. In Section 3 
we investigate new example of regularisation of a weak singular¬ 
ity. In Section 4 we discuss regularised by varying G cyclic and 
non-singular in matter density and pressure models of the universe. 
In Section 5 we study the multiverse model of constant entropy 
composed of cyclic individual universes of growing or diminish¬ 
ing entropy. In Section 6 we consider similar multiverse model but 
within the framework of Brans-Dicke theory of varying G replaced 
by a dynamical field 0. In Section 7 we give our conclusions. 


2 VARYING CONSTANTS COSMOLOGIES 


We consider Friedmann universes within the framework of vary¬ 
ing speed of light c theories and varying gravitational constant 
G theories. The field equations read as (Barrow 1999; Albrecht 
& Magueijo 1999; Barrow & Magueijo 1999; Gopakumar & Vi¬ 
jayagovindan 2001) 


d 
a 

where p is the pressure and g is the mass density, and the modified 
continuity equation is 
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Here a = a{t) is the scale factor, a dot means the derivative with 
respect to time t, G = G(t) is time-varying gravitational constant, 
c = c(t) is time-varying speed of light, and the curvature index 
k = 0, ±1. The idea to get Eqs. (2.1) and (2.2) is that one varies 
the generalized Einstein-Hilbert action in a preferred frame (since 
the Lorentz symmetry is violated), and then gets the Einstein field 
equations with c varying (playing the role of a scalar field) (Bar- 
row 1999; Albrecht & Magueijo 1999; Barrow & Magueijo 1999; 
Gopakumar & Vijayagovindan 2001) (but see also (Ellis & Uzan 
2005; Ellis 2007; Magueijo & Moffat 2008) ). The Eq. (2.3), on 


the other hand, can be obtained by direct differentiating (2.1) and 
inserting into (2.2). 

It is worth mentioning the similarity of varying c equations 
above with the field equations for Friedmann universes with vary¬ 
ing fine structure constant a (or alternatively, the electron charge 
e (Bekenstein 1982)) which are as follows (Barrow et al. 2004; 
Alexander et al. 2016) 
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and 


(2.4) 

(2.5) 


ifi + 3-ij/ = 0, 
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( 2 . 6 ) 


where g^ a stands for the density of radiation and 



stands for the density of the scalar field xj/ (standard with e = +1 and 
phantom with e = -1 (Caldwell 2002)) responsible for the variation 
of the fine structure constant 


a{t) = (2.8) 

It is easy to notice that (2.6) integrates to = C/a^. Also, one 
should mention that bearing in mind (2.5) and (2.7) one cannot get 
the universe with a bounce in the standard sense (a ’’non-singular 
bounce”) unless e = -1. In other words d < 0 throughout the time 
of the evolution. Putting this other way, for e = -t-1 there is no way 
to go from a contraction where H = a!a < 0 to an expansion where 
H = d/a > 0 through a singularity since H should have a minimum 
while in general one has 

(2.9) 

if the null energy condition is fulfilled. However, a ’’singular 
bounce” where there is no change of the sign of the second deriva¬ 
tive of the scale factor is still possible, and will be discusses in the 
next sections. 

Another point is that the theory represented by the Eqs. (2.1), 

(2.2) , and (2.3) is not the same as the Brans-Dicke theory in which 
G is replaced by a dynamical scalar field (f)(1) which possesses an 
extra field dynamics equation of motion (analogue of (2.6) in vary¬ 
ing a theory) and such a system is stronger constrained than (2.1)- 

(2.3) (see our Section 6). One of the points is that in both Brans- 
Dicke (Sakai & Barrow 2001; Bhattacharya et al. 2015) and G(l) 
theories (Mersini-Houghton 2014; Sidhart et al. 2015) based on 
(2.1)-(2.3) black holes other than those of general relativity cannot 
exist (”no-hair Brans-Dicke theorem”) though it is not necessarily 
the case in varying speed of light theories (Magueijo 2001). 


3 REGULARISING WEAK SINGULARITIES AND 
EACING A LITTLE-RIP 

In Djbrowski & Marosek (2013) we have shown how regularisation 
of the series of standard and exotic singularities (e.g. D^browski & 
Denkiewicz 2010) by varying constants was possible. In particular, 
it applied to sudden future singularities which due to regularisation 
could be gone through in a similar fashion as in D^browski (201 1). 
In fact, their weak nature was considered by Fernandez-Jambrina 
& Lazkoz (2004, 2006); Fernandez-Jambrina (2007, 2010). 

Here we use the scale factor which after appropriate choice of 
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parameters admits big-bang, big-rip, sudden future, finite scale fac¬ 
tor and w-singularities and reads as (Djbrowski & Marosek 2013) 


a(t) = flj 





(3.1) 


with the constants a,, m, n (e.g. Barrow et al. 1986; Barrow 2004; 
Freese et al. 2008; Shtanov & Sahni 2002), where we have split the 
scale factor into the two factors, one giving a big-bang singularity 
a SB, and another giving an exotic singularity We modify (3.1) 
into 
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ex 


Qs exp 
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which allows to extend the scale factor onto the times larger than tj 
(f > t.,) (Barrow & Cotsakis 2013): 


(3.3) 
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(3.7) 


for t > ts. There are various cases of weak singularities: for 1 < « < 
2 there is a sudden future singularity, for 0 < n < 1 there is a finite 
scale factor singularity, and for n > 2 there is a generalised sudden 
future singularity. 

Following Djbrowski & Marosek (2013) we may assume the 
variability of G as 


G(r) = Go 



(3.8) 


which allows the energy density (3.4), and the pressure (3.5) be 
regularized for r > 2 - n at r = f„ and the universe may con¬ 
tinue its evolution until t —» co, where the scale factor diverges 
(a —> oo) and both the mass density p ^ oo and pressure \p\ —> oo 
diverge, achieving a (strong) little-rip singularity (Frampton et al. 
2011) (which is a version of a (strong) big-rip (Caldwell et al. 2003) 
being reached in an infinite time). The regularisation of a sudden 
singularity due to (3.8) happens at the expense of getting an in¬ 
finitely strong gravitational coupling (G —> co) at f = ts, which is 
the characteristic being explored later in the paper. 

Another option to regularise a sudden future singularity 
(which admits a blow-up of pressure in (3.5) is to use the variahility 
of the speed of light 


c 


Co 


1-^ 
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(3.9) 


which happens for p > 2 - n at f = C though at the expense of light 


velocity slowing down to zero (c —» 0) at r = As in the previous 
case, after going through a regularised singularity, one faces a little 
rip for f —» oo, where a —» oo,p —> oo, and \p\ —> oo. Unfortunately, 
the little rip is a strong singularity and so we cannot continue the 
evolution of the universe further in order to get the cyclic evolution 
as in (Tolman 1934; Barrow & Djbrowski 1995). However, there is 
a transition through a sudden singularity as being a weak singular¬ 
ity which is regularised by varying constants (3.8) and (3.9). 

In summary, we have considered the scenario in which the 
universe starts with a constant scale factor with m = 0 in (3.1) 
(a bounce), goes through a sudden singularity, and then reaches 
another singularity - a little rip - found by Frampton et al. (2011). 
However, in this case we cannot continue the evolution towards the 
cyclic evolution since the little rip is a strong singularity. We will 
then try to find solutions which can bring the viable cyclic scenario. 


4 REGULARISING STRONG SINGULARITIES BY 
VARYING G 

As one can see from (2. l)-(2.2) and (2.4)-(2.5), the variability of the 
fundamental constants admits more freedom to a system of cosmo¬ 
logical equations. This eases a construction of the regularised cyclic 
models of the universe which will be shown below. 


4.1 sinusoidal pulse model - regularising big-bang 


In this section we use the equations (2.1) and (2.2) for positive cur¬ 
vature {k = - 1 - 1 ) and c = 0 in order to construct a fully regularised 
(in mass density and pressure) cyclic universe assuming that the 
scale factor takes the form 


a(t) = flfl 



(4.1) 


with Co =const., so the scale factor is "singular” (reaching zero 
as in the big-bang scenario). The main issue here is an additional 
assumption that the gravitational constant varies as (the speed of 
light c is taken to be constant here) 


G(r) = 


Go 

a\t)' 


(4.2) 


The plots of (4.1) and (4.2) are given in Fig. 1. It is easy to notice 
that the standard turning points with curvature singularities (zeros 
of the scale factor) are regularized by the infinitely strong gravita¬ 
tional coupling (G —> co). 

This can be seen by inserting (4.1) and (4.2) into (2.1) and 
(2.2) after which we obtain that the mass density p(t) and the pres¬ 
sure p{t) are nonsingular and oscillatory (see Fig. 2) 


P(0 

Pit) 
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(4.3) 

(4.4) 


Similar nonsingular in mass density and pressure cosmologies, 
though also with oscillating non-zero scale factor (a "non-singular 
bounce”), driven by the balance between the domain walls and 
the negative cosmological constant were obtained in Djbrowski 
(1996) and recently reconsidered as the "simple harmonic oscil¬ 
lator” (SHO) universe (Mithani & Vilenkin 2012, 2014; Graham 
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a(t) 



Figure 1. The cyclic scale factor (4.1) and the cyclic gravitational constant 
(4.2). 


et al. 2014). Introducing 

A 


Pa = 


dtrGo \ t\ 


+ c = pimtc) > 0, 


(4.5) 


where (m = 0,1,2,3...), one gets the equation of state for (4.3) and 
(4.4) in the form 


p(t) + p(t)c^ = Pa >0, 


(4.6) 


which tells us that the null energy condition is always fulfilled (a 
’’singular bounce” in the sense of (2.9)). The same is true for the 
weak energy condition which additionally requires (4.3) (p > 0). 
For the strong energy condition we have 

plus the condition (4.3) which means that this condition is fulfilled. 
In view of (2.2) it shows that the universe is decelerating {a < 0). In 
general, the strong energy condition means that gravity is attractive, 
but in our case its atractivity is overbalanced by the strong coupling 
of G(t) at the singularity (G —> co). Since 

3 / Ti^al 

P (™4) = = const., (4.8) 

8^Go \ tl j 

then it is easy to notice that at r = 0 the equation of state is of the 
network of cosmic strings type 

p(0) = -^-p(0)c\ (4.9) 


p(f) 




Figure 2. The nonsingular mass density (4.3) and pressure (4.4) for the 
models (4.1) and (4.2). 


At the turning points, we have: 


P (ntc) = 


3c2 


SttGo 

r2 

p(ntc) = -- 


: const.. 
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where n = 1 /2,3/2,5/2..., and the equation of state there is given 
by 

1 / 2n^al\ 

Pinta) ^--p(nta)\c^ -(4.12) 

In other words, the equation of state varies from the form of the 
cosmic strings (w = -1/3) given by (4.9) to the form (4.12) which 
depends on the parameters of the model uq, tc, and c (cf. Fig. 2). 
The scale factor is zero for t = mta and the big-bang singularity in 
each cycle is regularized (a —> 0, p —> const., p —> const.) at the 
expense of having strong gravitational coupling G —» co at r = mtc- 
Bearing in mind the form of the scale factor (4.1), we have 


cosTT-j-, mtc < t < (m + 1)4, 

costTp, (m -I- 1)4 < t < (m + 2)4, 

which means that it both may increase (a > 0) and decrease (d < 0). 
Same refers to the Hubble parameter 


an t 

H = - = — cotrr—, 

Cl tr tr 


(4.13) 


which is positive for m4 < t < (ra- 1 - 1)4 and negative for (m -l-1 )4 < 
t < (m + 2)4. The second derivative of the scale factor fulfills a 
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simple harmonic oscillator equation 

ii = ~a{t) (4.14) 

h 

SO that it is always negative (deceleration) and this is exactly what 
we see in the upper panel of Fig. 1. No turning point of evolution 
for the scale factor a(t) is possible (a "singular bounce” in the sense 
of the scale factor, but a "non-singular” bounce in the sense of the 
matter density g and pressure p). 

In fact, from (4.6) we can conclude that the null energy con¬ 
dition is fulfilled for the models under study. However, from the 
second Einstein equation (acceleration equation) we can easily see 
that 



and so there is no possibility for the universe to "bounce" i.e. to 
change from H = a!a < 0 (contraction) to H = aja > 0 (expan¬ 
sion) since H should have a minimum. In our model G is varying 
and this is why we have an extra term on the left-hand side of (4. 1 5) 
, but it does not help a proper "bounce” for the scale factor (geom¬ 
etry) since 


7T^ Cp- 

H = -—-- <0. 

Defining 


(4.16) 


G a 

Hg = -=-2- = -2H (4.17) 

G a 

we can see that Hq > 0 always accompanies H < 0 {a contrac¬ 
tion of space is always balanced by an expansion of gravity) which 
makes an apparently singular and sharp bounce regular in matter 
density and pressure due to the special form of running gravity. 

Similar situation appears in ekpyrotic/cyclic scenarios rooted 
in superstring/brane theories (Khoury et al. 2001, 2002; Steinhardt 
& Turok 2002; Khoury et al. 2004; Donagi et al. 2001; Turok & 
Steinhardt 2005), where one has some special coupling of gravity 
in the Lagrangian of a 4-dimensional theory (Einstein frame) as 
below 


S=jd^x^g 


- F(0) +A>P)(eR + 0.) 

lottG 2 


(4.18) 


where the potential has an explicit form as 

V(0) = 1 / 0 ( 1 - F{^\ F{^) oc , (4.19) 

and is the energy density of radiation while q„ is the energy 
density of matter, gj - dilaton/string coupling constant, (f) the scalar 
field, P(^) the scalar field coupling, and s, Vo are constants. 

This scenario is admissible only for the collision of the bound¬ 
ary branes and composes of infinitely many such collisions - each 
of them produces a big-bang after which the evolution repeats in 
contrast to a one brane collision in the earlier ekpyrotic scenario. 

The idea of a cyclic universe there is that the 5th dimension 
(orbifold) collapses, while the 4-dimensional theory has no singu¬ 
larity. This can be seen effectively by looking onto the generalised 
Friedmann equations which read as 

^ , (4.20) 

a = -'P+0'eR + \l^Q,i^ ■ (4.21) 
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The crucial point is that one chooses appropriate form of coupling 
jS oc 1/a, and we have 

1 1 


0R 


[o»)l‘ 


( 4 / 


const. , 


(4.22) 
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const. , 


(4.23) 


in the limit a —» 0, where the standard big-bang singularity usually 
appears. 

In our approach, the role of the coupling //(i^) oc 1 /a is played 
by the running gravitational constant G(t) oc l/a^(f) which regu¬ 
larises the mass density and pressure analogously as in (4.22) and 
(4.23) of cyclic scenarios. So one has analogously a kind of "sin¬ 
gular bounce” in the scale factor a(t) and a "non-singular” bounce 
in the mass density and pressure. 


4.2 tangential pulse model - regularising big-bang and 
big-rip 


Another interesting case is possible when one chooses the scale 
factor to be 


a{t) = Co 



/ ^ \ 

Co 

tan 7T— 

\ tsi 


(4.24) 


and the gravitational constant to vary as 


G(0 = 


4G, 

sin^ 


(4.25) 


We can see that the scale factor (4.24) is infinite (a —» oo) for 
t = nts with n = 1/2,3/2,5/2,... (like at big-rips) and it is zero 
for t = mts with m = 0,1,2,3,... (like at big-bangs) so that we can 
say that we face "singular bounces” (see Fig.3). Besides, each time 
the scale factor a(t) attains a singular value (vanishes or reaches in¬ 
finity), the gravitational coupling becomes infinite (G —> oo). This 
is somewhat strange since it regularises both big-bang and big-rip 
so its "attractive” or "repulsive” nature does not seem to have any 
difference here. 

The mass density and pressure are given by 


P(0 

P(t) 


3 

8;rG, 

•7 

C“ 

~87rG~, 


3c^ cos‘^ (n^) 

- 1 - —— 

al 

^ 


> 0, (4.26) 

C cos"^ 


(4.27) 


from which it is clear that they are regular (cf. Fig.4) - in that 
sense, we have regularised big-rips and big-bangs ("non-singular 
bounces”). From (4.26) and (4.27), one calculates that the equation 
of state is 


pip) = -C 




7T2 Co 
Ctp 


3G, 


4G/f2 


(4.28) 


where 

3k 

P > Pmm = pints) = ^77^. (4.29) 

which is in agreement with (4.26) for t = nt,. The maximal values 
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of the energy density and pressure are given by 


Qmax ~ QO^^s) ~ 


Pmax = P{mt,) = ■ 


8;rGj \ r? 


^ H-^ 


SttGs 


n c 

-r + -r 


(4.30) 

(4.31) 


The plot of (4.28) is given in Fig.5. 

Bearing in mind (4.26) and (4.27), one sees that the null en¬ 
ergy condition 


cV (t) + P it) = 


4;rG,, 


7T^ y 

-r -2- 


sin^ 


+ 4- 




(4.32) 


is satisfied at f = ml., and violated at r = nt,. We can then say that in 
this context (cf. the discussion around the formula (2.9)) we have 
’’singular bounces” at big-bangs and ’’non-singular bounces” at big¬ 
rips (null energy condition is violated by phantom). It proves that 
the former is a big-rip-like singularity and the latter is a big-bang- 
like singularity. The word ’’like” comes from the fact that the mass 
density and pressure are regular which is not the case at a big-bang 
and a big-rip. Perhaps better names would be ’’regularised” big-rip 
and ’’regularised” big-bang. It is interesting to note that the pressure 
at both of these situations p(mt^) and pints) is negative. 

The strong energy condition reads as 


cV(0 + 3p it) 


3c^ 

4^^ 





(4.33) 


so that 


c^p imts) + 3p imts) = 


c^p inis) + 3p ints) = 


3c^ 
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> 0 , 


< 0 . 


Further, bearing in mind that 




37tC^ 




(4.34) 

(4.35) 


(4.36) 


one concludes that the equation of state is that of phantom 


pints) = - ^Q(nts)c'^. (4.37) 

An interesting fact for the model with such an equation of 
state (though for the standard Einstein relativity) was found by 
Fernandez-Jambrina & Lazkoz (2004). In fact, it does not lead to 
the null geodesics incompleteness. In our case one has it effectively 
due to the influence of the special gravitational coupling of G(r) 
given by (4.25). It is also interesting to note that the model of sim¬ 
ilar equation of state violates the conditions of the Green and Wald 
theorem of averaging local inhomogeneities in the universe which 
predict the gravitational radiation contribution only (Green & Wald 
2011, 2013). 

Finally, from (4.37) one sees that the dominant energy condi¬ 
tion is violated as well for t = nts- This proves that one deals with 
phantom-like behaviour of the model at r = nt,. 

In this tangential pulse model again no turning point of the 
evolution for the scale factor ait) is possible (a ’’singular bounce” 
in the sense of the scale factor, but a ’’non-singular” bounce in the 
sense of the matter density g and pressure p). 


<i(<) 




Figure 3. The cyclic scale factor (4.24) and gravitational constant (4.25). 


5 CYCLIC MULTIVERSE 

In this section we start with the first law of thermodynamics (con¬ 
servation of energy) allowing the fundamental constants to vary and 
extend our discussion of regularised universes onto the model of the 
multiverse as a set of complementary universes evolving in a way 
which fulfils the condition of keeping the total entropy constant or 
growing i.e. obeying the second law of thermodynamics (entropy 
law). 


5.1 thermodynamics of varying c 


We first consider the varying c models (G = 0) starting from the 
first law of thermodynamics 


dE = TdS-pdV, (5.1) 

where E is the internal energy, T the temperature, S the entropy, p 
the pressure, and V the volume. Since 

E = mc^ = pVc^ (5.2) 

where p is the mass density and s = pc^ is the energy density, then 
we can take the differentials of E getting (Youm 2002) 

dE = c^pdV + Vdp -H 2pVcdc, (5.3) 


which inserted into (5.1) and taking the time derivative gives 

+ (5.4) 

On the other hand, the continuity equation for varying c models 
gives 



3kcc 

AnGa^ 


(5.5) 
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pw 




one can rewrite (5.6) assuming that the pressure and the mass den¬ 
sity are proportional 

p = wQC^ (5.9) 


as 


S =2 


pVc^ c 


(5.10) 


From (5.10) we note a very nice analogy to an ideal gas equation 
of state provided that we take 


pVC 

- = const. 


Nks 

w 


(5.11) 


where N is the number of particles, kg is the Boltzmann constant. 
After integrating (5.10), it is easy to define the entropy in a form 
which is analogous to the standard Gibbs relation (it is really in¬ 
teresting to note that in this case the volume of the phase space is 
proportional to the speed of light) 


Nks 

5(r) = 2^1n[Aoc(f)], 
w 


(5.12) 


where we have included an integration constant Aq which, without 
the loss of generality can be taken to be one. All this simplifies for 
flat k = 0 models since for them 'p = p, and using the barotropic 
equation of state p = wpC with the barotropic index w = const, we 
have 


2 pV 2 

S(t) = -^lnc(t) = -Nkslncd) . 
w 1 w 


(5.13) 


Figure 4. The nonsingular mass density (4.26) and pressure (4.27) for the 
models (4.24) and (4.25). 



Figure 5. The equation of state (4.28) for the model (4.24)-(4.25). The pres¬ 
sure is always negative and reaches the same initial value (4.29) for t = nt^ 
andn = \ 12,1,12,511, .... 


5.2 thermodynamics of varying G 


For the varying gravitational constant G (c = 0) models, one starts 
with the continuity equation (2.3) (which is an equivalent of the 
first law) to get 


V 

V 

so that 


P+y\P+-\ 


s = - 


pVC G 
~1~G' 



and again assuming an ideal gas equation of state 

pVC 

—— = const. = Nkn, 


one gets for varying G that 


(5.14) 


(5.15) 


(5.16) 


S(t) = Akflln 


Afl 

W) 


(5.17) 


where Aq is an integration constant (and can be taken to be one). 


Comparison of (5.4) and (5.5) leads to 


Vc2 

Defining 


S = ■ 


3kcc 

AnGa^ 


P = 


3 I cP kC 
^ttG \a^ 

1 kc* 


P = P~'5~F~2 
SttG 


= p + 


3 kC 


(5.6) 

(5.7) 

(5.8) 


5.3 varying c cyclic doubleverse 

From (5.12) and (5.17) one can see that the entropy is a function 
of the fundamental constants and it can increase or decrease de¬ 
pending on the form of the functions c(t) and G(f). According to 
the second law of thermodynamics entropy should never decrease 
and so we have to pick up some restricted form of these functions. 
Flowever, this should be the case in our universe only. In the mul- 
tiverse composed of individual universes this can be different. Our 
idea is that the total entropy of such an entity can still be constant 
or increase, while the entropies of individual universes can either 
increase or decrease. 
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Let us assume that the entropy of the multiverse is constant. 
Since the entropy is an extensive quantity, then we have 

n 

S =^5, = 5,+52 + 53 + ... + 5„=0. (5.18) 

1=1 

Specifically, for two universes which form a multiverse (which 
from now on we will call the “doubleverse”) from (5.13) we have: 

5, = ^7V,A:Bln[ci(f)], (5.19) 

w 

S2=lN2kBln[c2m (5.20) 

w 

We assume the following ansatze for the functions c{t): 

Cl it) = (5.21) 

C2 (r) = (5.22) 


where Ai and A 2 are constants. The total entropy of the doubleverse 
is 

Piy\c] P2V2C2 

S =Si+S2 = 2-^Ai<l,i+2^-^A2Mt) ■ (5.23) 

An interesting case which can be solved exactly is when 


Q\Vic\ 

Ti 





(5.24) 


or on the footing of (5.11) 


A/iAi — N2A2 , 


(5.25) 


c(I) 




which allows to pick up the following functional dependence of Cj 
and C 2 


(pi it) = sin^ j 

t)- 

i L/ 

(p2 it) = COS^ 

U) 

\ tj 


and write (5.19) and (5.20) as 


51 

52 


^ 2 
— N{kB^.\ sin 
w 



^ 2 

— N 2 kBA 2 cos 
w 



(5.26) 

(5.27) 


(5.28) 

(5.29) 


The plots of (5.21)-(5.22) and (5.28)-(5.29) are given in Fig. 6. 

An example of a doubleverse model which fulfils the field 
equations (5.7), (5.8), and (5.9) is the one with w = -1/3, 
which corresponds to already mentioned cosmic strings universe 
for which the scale factors are 


a{t) = ai{t) = fl 2(0 = Out. 


(5.30) 


On the other hand, assuming that still the pressure and the energy 
density are proportional 


Pit) 
ct it) 


~ pit). 


(5.31) 


we can extend (5.9) with vi> = -1/3 into the form 


Pit) 
ct it) 

where 


-3/5(0 + 


Poit) 
ct it) ’ 


Poit) ^ 1 a 

c^it) AttG a' 


(5.32) 


(5.33) 


Figure 6. The speeds of light (5.21), (5.22) and the entropies (5.28)-(5.29) 
of the varying c doubleverse. 


and so (5.1 1) transforms into 


ip - Po) F 
T 


= const. 


(5.34) 


which is an equation of state in the form of the van der Waals im¬ 
perfect gas model. 

Now, one may introduce the two scale factors equal and of the 
similar form as in (4.1): 


ait) = fl] (0 = 02 it) = Co 



(5.35) 


which satisfy (5.7), (5.8) and the equation of state (5.32) with a con¬ 
stant po/c^ = 7r/(4GtJ). Then, the evolution is evidently cyclic. The 
scale factor (5.30) is recovered from (5.33) provided that po = 0. 
The point here is that although the geometrical evolution of the 
’’parallel” universes is the same, this is not the case with the evo¬ 
lution of the physical constants Ci(t) and C2(0 which is evidently 
different. 


5.4 varying G cyclic doubleverse 


For varying G models we use (5.17) to get: 


c , r , 

51 + 5 2 - — Tfi — In 
r 1 


Gi (0 


P2V2C^ 

H-— In 


G2 (0 


(5.36) 


We choose the relation between the gravitational coupling constant 
G(f) and the scale factor af) (taken in the form given by (4.2)) in 
the first universe as 

G,(r) = GA/a]it), (5.37) 


MNRAS 000, 1-7? (2015) 


















Cyclic multiverses. 9 


0{t) 




Figure 7. Gravitational constants (5.37), (5.38) and the entropies (5.36) of 
the doubleverse for the scale factor (5.35). 


where a\{t) is the scale factor of the first universe, while in the 
second universe we choose that 

G2(1) = Gealdl (5.38) 

where 02 (f) is the scale factor of the second universe. In order to 
compare different evolutions of the varying physical constants Gj 
and G 2 , as the first attempt, we assume that both scale factors are 
equal i.e. 

aiit) - a 2 {t) = a{t). (5.39) 


In other words, both universes have the same geometrical evolution. 
In order to get the total entropy constant, we should again apply the 
condition similar to (5.24): 


PlVlC^— _ P2V2C^~ 

-I — - A') . 

r, T2 

where 




P^2 — ln/l2G 




(5.40) 


(5.41) 


Then, we have two universes with the total (net) entropy constant, 
though it is decreasing or increasing in each individual universe. 
Comparing the scale factors in both universes we have that 

G,(i)^-^ (5.42) 

02(0 

and 

Si(f) oc lnG|(f) « In —oc(5.43) 
G 2 (t) 02(0 

This means that we have the same geometrical evolution of 


the two “parallel” universes, but their physical evolution (physi¬ 
cal coupling constants) is different. More precisely, the first uni¬ 
verse and the second universe have the same scale factors but their 
gravitational constants are different. In one of these universes the 
evolution of the gravitational constant is singular depending on the 
evolution of the scale factor a{t). This is an alternative approach 
than in the case of the Kaluza-Klein theory in which microscopic 
dimensions shrink and have the scale factor which diminishes while 
macroscopic dimensions expand (Polchinski 1998). 


5.5 discussion 

For the whole body of this Section it is also worth mentioning that 
the scale factors for each universe within the multiverse satisfy the 
field equations (2.1) and (2.2). 

From Fig.7, which has been plotted under the assumption that 
Ga = Gb and Qq = 1. one can see that the universe 1 may replace 
its evolution along the trajectory of the universe 2 at the maximum 
expansion point. It is an interesting new option for the evolution 
of the universe - now put in the context of the multiverse - where 
some effects take place at the turning point. Similar effects (though 
considered to be appropriate to the same universe) were studied 
in the context of quantum cosmology and the probability of tun¬ 
neling from an expanding phase into a harmonic oscillator phase 
was calculated (D^browski & Larsen 1995). This discussion was 
more recently developed under the name of the simple harmonic 
universe (SHU) scenario (Mithani & Vilenkin 2012, 2014). All this 
is in agreement with the claim that macroscopic quantum effects 
in cosmology are possible (Gold 1962; Kiefer 1988; Kiefer & Zeh 
1995; D^browski et al. 2006). 

The same is possible for the varying c multiverse of subsection 
5.3 where at the maximum expansion of the scale factors (5.35) the 
universes may replace their evolution due to quantum effects. 

Here we mean the multiverse which is quantum mechanically 
entangled and there are periods of its evolution where this entan¬ 
glement matters (e.g. the maximum expansion point) turning and 
enters the behaviour of individual universes, but most of the evolu¬ 
tion of the individual universes is classical. In other words, we try 
to say something about possible parallel evolutions of the two clas¬ 
sical patches of the multiverse (which can be understood on some 
lower levels of the hierarchy given by Tegmark (2003) - from I to 
IV). These two patches are classically disconnected, but they can be 
quantum mechanically entangled (Robles-Perez & Gonzalez-DIaz 
2010; Robles-Perez & Gonzalez-Diaz 2014; Robles-Perez et al. 
2015) and the effect of entanglement can be imprinted in one of 
the universes (for example in CMB) or become explicit when these 
universes reach some special conditions (as we have proposed for 
example - at the maximum expansion). 


6 CYCLIC BRANS-DICKE MULTIVERSE 

Brans-Dicke field equations (Brans & Dicke 1961) can be written 
down as (and following Section 5 are valid in each individual uni¬ 
verse of the multiverse separately giving the same form of the field 
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equations ( 2 . 1 ) and ( 2 . 2 ) in the limit (j) = const.) 


pit) = 


30(f)[a2(r) + fcc2] 3a (f) 0(f) atifr it) 


Pit) = 


it) Stta it) 

(p it) \a^ it) + kc^ 


16;r0(f) 

2a it) 0 it) (ji>(j)t (t) 


( 6 . 1 ) 


□0 = 


2<p it) a it) 
a (t) 

8tt 


(t) 


- 0 (f) 


a it) 


Imp (f) 


3 + 2£j 


( 6 . 2 ) 

(6.3) 


and the last of them (6.3), absent in the system (2.1)-(2.3), is equiv¬ 
alent to the continuity equation in Einstein relativity 




(6.4) 


entropy. It also differs from the conservation equation of varying 
G theory (2.3) since it does not contain the term —qGIG which 
couples mass density with variation of gravitational interaction. 

However, if we want to assume that the entropy is not constant, 
then we need an extra entropic term in the continuity equation, get¬ 
ting 


. Vl P\ ,.TS<p ^ 


(6.5) 


This leads to a modification of the field equation (6.3) in Brans- 
Dicke theory by an additional function (we take k = -l-l in (6.1)- 
( 6 . 2 )) 


= 


8;r 

3 -I- 2oj 


T + Nit), 


( 6 . 6 ) 


and the function N it) is responsible for an increase or a decrease 
of the entropy. 

Such modification has recently been considered by Kofinas 
(2015) in the context of energy exchange between a scalar field and 
the standard matter. This allows an avoidance of having the prob¬ 
lem with singular conformal cosmology (and also cyclic conformal 
cosmology - CCC (Penrose 2009, 2010; Majid et al. 2008)) limit 
o) = -3/2 (or Kofinas’ A = 2/(3-H2ai) —> oo) and allows the entropy 
growth needed to drive our cyclic multiverse models (Blaschke & 
D^browski 2012). 

Comparing (6.5) and ( 6 . 6 ), one gets 


S = 


Nit)V P 
IbnT p 


(6.7) 


Since the mass density and the pressure are both the functions of 
time, then we can make N it) the function the of the mass density 
and the pressure N(p,p) and assume that 


Nip,p)V 

T 


■ const. 


( 6 . 8 ) 


For variable entropy, one modifies the continuity equation (6.4) to 
the following form 


.. 3d (f) 0 (f) 8;r [p it) - 3p (f)] 

p H-7^^- = -;;-;;-h A'(f), 


aif) 


3 -t- 2u) 


where 


M(f\ - it) ^ttfjp jt) 

3 -I- 2£j 3 -I- 2m ’ 

and a and p are some constants. 


(6.9) 


( 6 . 10 ) 


In such a case the entropy equation (6.7) integrates to give 

( 6 . 11 ) 




IT (3 -H 2m) 

where for a 0 we need to use an equation of state p ~ p. 

In a special case of the low-energy-effective-string theory 
(Gasperini & Veneziano 2003; Lidsey et al. 2000), where m = -1 
and k = -l-l, one can take the ansatz (4.1) in the first universe (of 
the doubleverse) as 


( 6 . 12 ) 


fll (t) = Oa 

■ 1 t\ 

sm 7T — 


\ tj 


and 


01 (0 ~ a\ it) = Pa sin | 7 r— 
ts 


(6.13) 


where = 4/5, Pi = -32/15 and the speed of light c = { 0 ^ 71 )!itP 
necessarily in order to fulfill the field equations ( 6 .l)-(6.3). In this 
case the mass density and the pressure read as 




3 -I- 11 cos^ (;r — 


piit) = 


C^pAtt 

StP 


■ 9 cos \7T — 


while the equation of state is 

Piit) _ 9 150A?r 

—--Yyp> « + 


(6.14) 

(6.15) 

(6.16) 


In the second universe (out of the doubleverse), we can assume that 


flz it) = Ob 


and 


p2 it) ■ 


. 

t\ 

sin 

h) 


til it) sin^(;r^) 


(6.17) 


(6.18) 


with 012 = 0 and P 2 = -4 and the speed of light c = (agtrj/f,. In this 
model, the mass density and the pressure take the form 




2 -I- sin^ {.ttf) 


P2 it) 


C^p27t 

■”8C^ 


sin"* ijt-f) 

2 - sin^ [ttp] 


sin"* (?tp) 


(6.19) 


( 6 . 20 ) 


and we used the fact that for 02 = 0 the equation of state is not 
necessary. Finally, provided that 

1 j. [®iPi (^)], . [®2P2 ( 0 +//2P2 ( 0 ] V2 

)ttpA --=ln^B- - -, (6.21) 

11 I 2 

the sum of entropies of individual universes in such a doubleverse 
is constant i.e. 5=0 (cf. Eq. 5.18), though it changes in each of 
them. 

As we can see from (6.12), (6.17) and (6.13, (6.18) the ge¬ 
ometrical evolution of the scale factors is similar and cyclic as in 
the previous models of Sections 4 and 5, while the evolution of the 
physical constants represented here by the dynamical fields pi and 
p 2 which play the role of Gj and G 2 can be different. In other words, 
all the discussion of the physical effects of Section 5 is valid. 
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7 CONCLUSIONS 

In this paper we have extended our idea of regularisation of cosmo¬ 
logical singularities due to the variability of the fundamental con¬ 
stants (D^browski & Marosek 2013) onto the cyclic universe and 
multiverse models. We have shown new examples of regularisation 
of the sudden future singularity by varying c and G though continu¬ 
ing towards a little-rip singularity which prevented the construction 
of cyclic universes. 

Then, by an appropriate choice of the evolution of the scale 
factor a{t) and the dynamics of G(t), we have been able to construct 
cyclic models in which the energy density and the pressure are non¬ 
singular and oscillatory. 

Being encouraged by the success of such a scenario we have 
extended this idea onto the multiverse containing cyclic individual 
universes with either growing or decreasing entropy, though leav¬ 
ing the net entropy constant. In particular, we have considered the 
“doubleverse” with the same geometrical evolution of the two “par¬ 
allel” universes with the evolution of their physical coupling con¬ 
stants c(l) and G(t) being different. We have found an interesting 
possibility for one universe (of the doubleverse) to exchange its 
evolution history with another universe at the turning point of ex¬ 
pansion. Similar investigations have been made in the context of 
quantum cosmology, where the universe (considered as the same 
object/entity) could quantum mechanically tunnel into an oscilla¬ 
tory (harmonic) regime. 

We have also found that similar scenario is possible within the 
framework of Brans-Dicke theory of varying G(l) oc 1 /0(t), where 

is the Brans-Dicke field though these two theories differ in the 
form of the continuity equation. 

In fact, we face various categories of the "bounces”. We have 
differentiated the "singular bounces” and "non-singular-bounces” 
referring to the appropriate quantities such as the scale factor, mass 
density and pressure. We have also calculated the energy condi¬ 
tions and noticed that the null energy condition is not necessarily 
violated due to the regularisation by variability of G. The discus¬ 
sion of the energy conditions shows that the null energy condition is 
fulfilled for the big-bang type bounces and violated for the big-rip 
type bounces. 

Our discussion is heuristic and preliminary so we are not able 
to address all the physical problems which are related to the issue 
of cyclic models within this scenario. Besides, we stick to purely 
isotropic geometry models not discussing the problem of possible 
anisotropy domination which can be solved by the very stiff equa- 
fion of state (p = wp, w » 1) during the collapse phase. We have 
explicitly found analogies with ekpyrotic/cyclic approach where 
the role of varying G is played by the appropriate coupling of mat¬ 
ter with a scalar field. The ghost (phantom) instabilities may not be 
faced here (at least partially) mainly because in most of the cases 
still there is no null energy condition violation. 

Finally, it should be said that the main problem in the liter¬ 
ature related to the multiverse idea concentrates on the measure 
problem, but we have not investigated it in any aspect here (see 
e.g. Lehners (2012)). Instead, we try to say something about pos¬ 
sible parallel evolutions of the two classical patches of the multi¬ 
verse (which can be understood on some lower levels of the hier¬ 
archy I to IV of Tegmark (2003)). Individual universes/patches of 
the universe are classically disconnected, but they can be quantum 
mechanically entangled and the effect of entanglement can in prin¬ 
ciple be imprinted in individual universes or become explicit when 
these universes reach some special conditions. 
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